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ABSTRACT 

Using the theoretical formalism of Aim and Schuck, we have 
studied the behavior of collective modes over the whole coupling 
regime from weak coupling (high density) to strong coupling (low 
density). We observe that in the weak-coupling limit particle-hole 
RPA is approached to low momenta and the collective mode in the 
strong coupling limit reproduces the Bogoliubov mode for the 
weakly interacting bosons. 
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INTRODUCTION 

The one dimensional Fermi gas with 
attractive 5 interaction among the fermions 
has been used as a model to know the 
physical properties of realistic Fermi system. 
The exact solution for its ground state 
energy is known from the Bathe ansatz 1 , 
therefore one can test approximate solutions 
of the systems. Quick, Esbagg and deLalno 2 
have developed mean field approximation 
basically of two types (a) plane wave Hartre- 
Fock (HF) (b) non-plane wave (HF) and 
BCS. It was found that the BCS solution can 
describe the cross-over between weak 
coupling (as a weakly interacting gas of 
fermions) and strong coupling (gas of 
bosonic two particle pairs) in the system. In 
particular, they find that the BCS solution 
for the ground state energy coincides with 
the exact solution in both weak and strong 
coupling. In this sense, the system may serve 
as a simple model to study the transition 



between weak and strong coupling 
superconductivity in ID Fermi system. This 
transition between weak and strong coupling 
has been discovered by Leggett 3 and by 
Noziers and Schmitt Rink 4 in three 
dimensions. The simple form of the 
interaction allows one to carry out 
approximation beyond the mean-field level 
such as ordinary random phase 
approximation (RPA's) or generalized RPA's 
in controlled fashion. In this approximation, 
one is able to calculate contribution to the 
ground state energy of the system beyond 
the mean field ground state energy. Williams 
and Bloch 5 have also discussed the ordinary 
(particle hole) RPA for the one dimensional 
electron gas. The same approach was also 
applied by Bernner and Haug 6 . Friesen and 
Bergersen 7 have applied Singwi Sjolander 
generalization of RPA 8 to the ID electron 
gas. This model can serve as an initial 
approximation to realistic systems such as 
quasi one dimensional metal 9 . 
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In this paper, using the theoretical 
formalism of T. Aim and P. Schuck 10 , we 
have studied the collective modes of the one 
dimensional Fermi gas within the quasi 
particle random phase approximation (RPA). 
One calculates the collective excitations for 
the ID attractive Fermi gas with a 8 
interaction at T = 0 by applying the quasi 
particle RPA. One obtains quasi particle 
RPA equation for the two quasi particle 
propagator starting from the Bogoliubov 
transformed Hamiltonian of the system 
using the equation of motion for the two 
particle Green's function. The homogeneous 
two particle equation yields the condition for 
the collective excitation in the system. It 
coincides with result found by Anderson 11 , 
Rickayzen 12 , Bardasis and Schrieffer 13 and 
others 14 ' 15 by the equation of motion method. 

MATHEMATICAL FORMULAE USED 
IN THE STUDY 

One writes the Hamiltonian for one 
dimensional Fermi gas in second 
quantization 

H = 2> + T I <*.*#! K 3 K 4 )C ki C+ Q 3 Q 4 

(1) 

where k denote momentum and spin 
quantum number of the particle and 

(kfa \V\K 3 K 4 ) is the ant i_ symmetr i ze d 

matrix element of the two body interaction. 
Using Bogoliubov transformation of the 
creation and annihilation operators 
Hamiltonian (1) is transformed into a new 
form. The transformed Hamiltonian is 
derived by several author, 1618 The 
transformed Hamiltonian is written in the 
following form 



k\k^ ^"1^2 

(2) 

In equation (2), H° is the BCS 
ground state energy and H 11 is the diagonal 
part, H 20 is the off-diagonal part of the 
Hamiltonian. If one demand that the off- 
diagonal part H 20 vanishes then one obtains 
the well known relations of BCS gap 
equation 19 . 

The H 31 term in the Hamiltonian in 
eq n (2) does not contribute to the RPA 
equations. The other terms H 40 and H 22 in 
the Hamiltonian describe the residual 
interaction among the quasi particle. 

These terms are neglected in the 
model BCS approximation. Therefore in 
order to go beyond the BCS mean field 
approximation one has to include the 
residual interaction among the quasi particle. 

Now, one can treat the quasi particle 
within the generalized RPA approximation. 
For this, one introduces two particles 
Green's function with respect to quasi 
particle basis and derive the equation of 
motion for the Green's function G. The 
equation has the form of Dyson equation 20 
for the two particle propagator matrix G. 
The static elements of the mass operator A 
and B are given by the double 
commutation 21 
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K l K 2 K 3 K 4 KjK 2 K 3 K 4 \ K l K 2 / K 1 K 2 K 3 K 4 ,-,11 / , , / \ ^21 / 1 1 f \ 

G [kkq,G))—G \kkq,G)) 



\ W4 =(K^ 2 ,[^,«]] 



and 



B 



k] k<jk-ik,A 



(3) 



(4) 



= Gl 1 (kq, 0))S kX + , . 
+fin(fc, 9 )A 4 +£ ii /(fc^)r tv 



[2flw(*, ? )]Z tv 



(7) 



0 



The quantities Z, A and r are given by the 



Here \ ' refers to averaging with respect to 
the BCS ground state. But in place of BCS ec l uatlons 
ground state one has correlated ground state 

which corresponds to a generalized quasi Z r?ffl = --^m(A:^)[G 11 (r%«)+G 21 (r%«)] 



particle RPA . 

Taking 5 interaction one solves the 



(8a) 



system equation in the CO representation. One A* q0J =-vY]{k'\q\G l \k'T<q,0^+G 2 \k'T<q,0^ 
obtains coupling equation for G 11 and G 21 ' ' '* 1 



k" 



also for G and G . One obtains 



(8b) 



Gl\kq,0)) = 
Gl 2 (kq,0)) = 



1 



k,q,a> 



= -vYj{k\q\G l \k'M,(o)+G 2 \k''kq,0^\ 



(8c) 

The quantities m(k,q), n(k,q) and l(k,q) are 
combination of u k and v k and is given by 



(5) 



where 



Ek, q - E k + E k+q 



We have assumed that collective 
pairs have zero total spin One arrives on the 
expressions 



m(k,q) = u k v k+q + v k u k+q 

n(k,q) = u k u k+q +v k v k+q 

£(k,q) = u k u k+q -v k v k+q 
and 

1 



G n (kk'q,0)) + G 1 \kk'q,o)) 
= Gl\kq,ai)S Ki ,+- 



+E k Xk,q)A +a*{k,q)r„ q „ 



[ 2E k, q n < k ><l)] Z k' q ,a 



Where 



1 + ^- 



(9a) 
(9b) 
(9c) 

(9d) 



HF 



£ k =^k+v k -M 



(6) 



(9e) 
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The quantity I a b c are in the notation 



Multiplying in equation (7) 
subsequently with m with n and in equation 

(7) with ^ and summing over k, one arrives 
at the system of equation for the quantities 
Z, A and r 

vI m A<lM 1+v/ w(?' <y ) 2v U(«- ffl ) 



-vn(k,q)Gl l 
-vi{k,q)Gl l 
-v/2m(k,q)Gl x 



x 



7 



^a,b,c ^ j 



a(k,q)b{k,q)c{k,q) 



en) 



(10) 



Here a(k,q) = [E kq co] and b(k,q), c(k,q) = 

[n(k,q), ^ (k,q), m(k,q)]. This is a linear in 
homogenous system of equation for the 
quantities A, T and Z. This can be solved 
with the help of matrix inversion. Equation 
(10) is a eigen value problem for the 
determination of the collective modes in the 
quasi particle RPA. The condition for the 
nontrivial solution is the vanishing of 
determinant. 



l + v/ £ ,„,„ {q,ti) vl^ nj {q,a) vI E ^ m {q,£l) 

vIn,„A<l>n) l + vI EJJ {q,£l) 2vI nfm {q,Q) = 0 
v/2I E , n , m (q,&) v/lI^iq.O) \ + vI Em J q ,Q) 
where Q. (q) denotes the eigen value for the collective excitation. 



(12) 



Weak Coupling Case 

The analysis of the weak coupling 
case has been given by Belkker and 
Randeria 23 . One quotes their result for the 
weak coupling collective made in one 
dimension. 



Q(q) = 



l-v 



1 m 



-,V2 



m 



(13) 



where the density of states in one dimension 
for parabolic dispersion was used. The long 
wavelength collective modes in the weak 



coupling case have a phonon like spectrum 
and are independent of the gap. If one 
compares equation (13) with the small-q 
expansion of the collective modes in the 
particle hole RPA, one finds that both 
coincides. This means that in weak coupling 
the behavior of collective modes for small q 
is not changed from the normal particle hole 
RPA. 

However, for large q and in 
particular near the point q=2k F the quasi- 
particle differs from the particle hole RPA. 

Strong Coupling Limit 

Now, from equation (12) one can 
expand the determinant of the collective 
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modes for small -q and - Q. limit with 
respect to small gap. In the strong coupling 
limit the gap goes to zero with density. 2 One 
expands the gap equation in terms of A 2 



1 



1 + 



V 



2 

1 I J 



2tc{ \£ t 



[dk 



2n 



(14) 



dimensionality of the system. An analogous 
expansion of the BCS density equation 
yields 

n = — A 2 J 2 

In ~ (16 ) 

Now, the effective chemical potential can be 
expressed in terms of the density and the 
coupling strength as 



(//*) 



V2 



4i 



-m 



V2 v±, 



— mv — nv 
64 8 



(17) 



Here one has introduced integral Ji. It has 
been evaluated for the strong coupling limit 
i.e. for |i *<0. Here |i * is effective chemical 
potential including the quasi particle shift 

j. =jdk=- 

0 



f k 2 



2m 



+ JU 



J i2 =\dk 



2m 



Then 



+ jU 



1/2 



J l = mn 1/2 /[2\ju' 
J 2 =m ll2 7tl(2\jU^ 12 
J 3 ^3m l ' 2 ^/(2 1 ' 2 \ J uf 2 



(15a) 



(15b) 



(15c) 



Now in equation (14) the integral in 
the first line is convergent for a contact 
interaction without a cutoff due to the one 



In the limit of zero density or zero gap, 
respectively equation (17) yields the 
condition 

.2 



2\ju : 



mv 



(18) 



This shows that in the extreme strong 
coupling limit the chemical potential i.e. the 
energy to remove a particle from the system 
in just half the two particle binding energy - 
E 0 in the Vacuum. 2 

The long wavelength dispersion relation in 
the strong coupling limit is given by 



a(q) = cq 



(19a) 



where the sound velocity c depends on the 
gap A . Now using the low density 
expansion of the BCS density (16) to 
substitute A 2 by the density n, one obtain the 
expression for the collective modes in the 
strong coupling limit. 



a{q) = 



vn 



Am 



1/2 



(19b) 
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Introducing the pair mass m B = 2m and the 
pair density n B = n/2 the above expression 
reduces to 



4m 



1/2 



(19c) 



Equation 19(c) is the well known 
Bogoliubov dispersion relation for the 
weakly interacting Bose gas 24 ' 25 in the limit 



of small q which is linear in q i.e. phonon 
like. Thus starting from interacting fermions 
with an attractive interaction the quasi 
particle RPA in the strong coupling limit 
yields the dispersion relation for weakly 
interesting gas of bosons (two particle bound 
states). The magnitude of the repulsive 
interaction among the bosons in equation 
19(c) is given by the fermionic interaction 
strength v and is consistent with the result of 
Hanssmann 26 and others. 27 32 



Table T x 

An Evauated results of collective mode Q. for small momenta q in the weak coupling limit 
(|i*=3.386 E 0 ) The numerical solution is compared with weak coupling expansion (13). 



q/k F 


QJ\i 


Numerical solution 
Equation (12) 


Weak coupling 
expansion eq. (13) 


0 


0.00 


0.00 


0.1 


0.126 


0.132 


0.2 


0.268 


0.275 


0.3 


0.355 


0.368 


0.4 


0.568 


0.599 


0.5 


0.573 


0.702 


0.6 


0.626 


0.925 


0.65 


0.645 


1.027 


0.70 


0.689 


1.146 


0.75 


0.735 


1.207 


0.80 


0.794 


1.258 
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Table T 2 

An Evauated results of collective mode Q. for small momenta q in the strong coupling limit 
(|X*=(-)0.4966 E 0 ) The numerical solution[eq.(12)] is compared to the strong coupling expansion 

[eq.(19c)] and free particle dispersion 
£2 =q 2 /2m B (a 0 = Bohr radius) 





£2/E 0 


q/a 0 


Numerical 
Solution 

Eq.(12) 


Strong 
expansion 
eq.(19c) 


Free particle 
dispersion 
£2 = q 2 /2m B 


o 


0.000 


0.000 


0.000 


U.Uz 


0.0022 


0.0020 


0.0002 


0.04 


0.0032 


0.0032 


0.0004 


0.05 


0.0034 


0.0035 


0.0005 


0.06 


0.0037 


0.00355 


0.00282 


0.07 


0.0038 


0.00362 


0.00294 


0.08 


0.0042 


0.00375 


0.00307 


0.09 


0.0048 


0.00382 


0.00322 


0.10 


0.0052 


0.00397 


0.00375 


0.11 


0.0056 


0.00416 


0.00400 


0.12 


0.0062 


0.00438 


0.00425 


0.15 


0.0094 


0.00469 


0.00473 


0.20 


0.0126 


0.00525 


0.00098 



DISCUSSION OF RESULTS 

In this paper, using the theoretical 
formalism of T. Aim and P. Schuck, 10 we 
have studied the behavior of collective 
modes over the coupling regime weak 
coupling (high density) to strong coupling 
(low density). We observe that the treatment 
of the residual interaction in the Hamiltonian 
(2) within the quasi particle or generalized 



RPA allows one to study the behavior of the 
collective over the whole coupling range. 
Collective modes in one dimension were 
found for the case of an attractive 8 
interaction. One observes that in the weak 
coupling limit one recovers Anderson's 
results 11 whereas in the strong coupling limit 
the Bogoliubov dispersion relation 24 for the 
interacting Bose gas of two particle pairs can 
be derived from the quasi particle RPA. This 
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is consistent with the fact that the BCS 
theory is capable of describing the extreme 
strong coupling limit i.e. the gas of two 
particle bound state properly and reproduces 
the exact result for the ground state energy 
in this limit. 2 In table T i' we have presented 
for small momenta q in the weak coupling 
limit (|J*=3.386E 0 ). The numerical solution 
equation (12) is compared with weak 
coupling expansion equation (13). Both 
solution coincides at q/kF = 0.4. This shows 
the consistency of the numerical solution 
with the well known weak coupling result 
which was obtained by Anderson 11 in the 3D 
case. In table T 2 , we have presented the 
evaluated results of collective mode Q. for 
small momentum q in the strong coupling 
case = (-)0.4966 E 0 )' The numerical 
solution equation (12) is compared with the 
strong coupling expansion (equation 19c) 
and also with free particle dispersion 
£2 =q 2 /2m B . We observe that the numerical 
solution starts linearly in q and is consistent 
with the strong coupling expansion gives in 
equation 19(c). This confirms the 
interpretation of the collective excitations in 
the strong coupling limit as Bogoliubov 
sound modes of the two particle Bose gas 
that is formed in the limit. The free particle 
solution reached the full solution for large q. 
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